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Abstract 

We present a Dalitz plot analysis of the decay — > K^t^^t:^ based on 281 pb~^ of e^e~ 
collision data produced at the V(3770) by CESR and observed with the CLEO-c detector. We 
select 67086 candidate events with a small, ~1.1%, background for this analysis. When using 
a simple isobar model our results are consistent with the previous measurements done by E791. 
Since our sample is considerably larger wc can explore alternative models. We find better agreement 
with data when we include an isospin-two vr^vr'^ S-wave contribution. We apply a quasi model- 
independent partial wave analysis and measure the amplitude and phase of the Ktt and tt+tt"'" S 
waves in the range of invariant masses from the threshold to the maximum in this decay. 



*Submitted to The 2007 Europhysics Conference on High Energy Physics, July 19-25, 2007, Manchester, 
England. 
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In comparison to many other meson modes, the K~ti^ti^ decay is unique in 

many aspects. A large contribution of over 60% from a Kt^ S-wave intermediate state has 
been observed in earher experiments [H, 0, [s], 0, [sl, [6| . The large branching fraction, B{D~^ 
K^TT^TT^) = (9.51 ± 0.34)%, for this mode makes it the usual choice for normalization of 
other meson decay rates. Understanding its peculiar intermediate substructure will be 
very beneficial. Two identical pions in this decay should obey Bose symmetry. Assuming 
quasi two-body nature of this decay, it is dominated by the two identical sets of K~7r+ waves 
interfering with each other. This symmetry significantly reduces the degrees of freedom in 
the regular Dalitz plot analysis and allows us to apply a model-independent partial wave 
analysis. One would also expect a small contribution from the isospin-two 7r+7r+ S wave, 
which can also be explored in this decay. 

The K~i[^T[^ decay has been previously studied with the Dalitz plot technique in 

many experiments, including MARK III NA14 j2|, E691 Q, E687 ||, and E791 gQ. To 
achieve good agreement with their data, E791 \5\ has included a low-mass K~7i^ resonance, 
K, that significantly re-distributed all fit fractions observed in earlier experiments. Due to 
the difference in model, the results of this high-statistics experiment have been excluded 
from the world average by the PDG. Compared to previous experiments CLEO-c 0] has 
many advantages. Our sample is four times larger than E791, and it is very clean; the 
background is about 1% of the selected sample. The invariant mass resolution in this three- 
track D-meson decay is very good, ~2.5 MeV/c^. Our experimental conditions are similar 
to that of MARK III, where D mesons are produced with small momentum. The other 
experiments, all fixed-target, study the D mesons produced with large momentum, typically 
tens of GeV/c. 

Interest in this decay was outlined in Refs. [sl, 0, 13], which focus on a low-mass Ktt 
S-wave contribution. These authors reanalyze the E791 data with their own models. E791 
reinterpreted their own data with a model independent partial wave analysis and it is 
this approach we apply in our analysis with minor modifications. 

We use 281 pb~^ of e~^e~ collisions at y/s ^ 3770 MeV, produced by the Cornell Electron 
Storage Ring (CESR) and accumulated by the CLEO-c detector. This sample corresponds 
to the production of ~0.8 M D^D' pairs in the process e^e~ ^(3770) D^D~ . We 
reconstruct a one side D-meson decay using three tracks and applying the standard CLEO-c 



algorithms similar to Ref. [ll|. In order to select K 7r+7r+ decays, we use two signal 



variables, /S.E = Eo — Ebeam, and m^c = y -^beam ~ where E^eam is the beam energy, and 
E]:, and are the energy and momentum of the reconstructed D+-meson candidate. We 
require |A_E| < 2a{^E) and \mBc — ^d\ < '2'(j{mBc), where resolutions a{AE) = 6 MeV 
[(5.74±0.02) MeV in single Gaussian fit], and a{mBc) = 1-5 MeV/c^ [(1.36±0.01) MeV/c^ 
in single Gaussian fit] represent the widths of the signal peak in the 2D distribution and 
projections shown in Fig. [TJ We select 67086 candidate events for the Dalitz plot analysis 
1^. The fraction of background, ~1.1%, in this sample is estimated from the fit to rriBc 



spectrum of Fig. [H 

In this analysis we apply a formalism similar to that used in other CLEO Dalitz plot 



analyses [13|, [IJ] . The Dalitz plot is symmetric under the interchange of like-sign pions, so 
we do the analysis in the two dimensions of high Kn mass squared, m'^{K7T)high, versus low 
Kti mass squared, m'^{K7i)iow This choice folds all data onto the top part of the Dalitz 
plot, as shown in Fig. [2](left). Figure [2] also shows the m'^{Kn) (two entries per event) and 
m^(7r7r) projections. 

The efficiency across the Dalitz plot is modeled with simulated events that are fit to 
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Data: D+ K" 7r+7t+ 



Data: D+ K" 71+71+ 



3540707-031 



Data: D+ ^ K" 7r+7t+ 



o 

^ 1.87 


CD 

"1.86 



1.84- 
-0, 







ENTRIES 130799 


































































































































1 -0.05 0.05 0.1 

AE (GeV) 



-0.1 -0.05 0.05 0.1 

AE (GeV) 




1.86 1.1 
rriBc (GeV/c^) 



FIG. 1: Distribution of niBc ''^s. AE (left) for K^ii^ii^ signal candidates in CLEO data. 

Events from the shown ±2(7 bands crossing the signal region are plotted as /S.E (middle) and m^c 
(right) distributions. The background events are selected from the sideband hatched box, which 
has about the same range of K~ti^'k^ invariant mass as the signal box. 
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FIG. 2: Dalitz plot for data and its projections on m?(K 7r"'")(two entries per event) and m?{TT~^7r~^) 
axes. 



a two-dimensional third-order polynomial. There is a notable decrease of the efficiency in 
the corners of the Dalitz plot, which is modeled by the multiplicative sine- wave threshold 
functions. The background shape is estimated from the ttibc and AE sidebands. Many 
possible resonances can contribute to the decay, and a total of 7 different contributions are 
considered. Parameters describing these resonances are taken from previous experiments [l^ 
or phenomenological publications. Only contributions with an amplitude significant at more 
than three standard deviations are said to be observed. Contributions that are not significant 
are not included in the decay model used for the result. We estimate a fit quality using 
Pearson's calculated for adaptive bins on the Dalitz plot. 

First we compare our results, obtained in the framework of isobar model (IM), with E791 
Models A, B, and C from Ref. [5|. In particular, the most advanced Model C contains 
K*{892)n+, i?*(1430)7r+, K*{1680)7i+, i^*(1430)7r+, i?*(1680)7r+, k7i+, and non-resonant 
contributions. Following E791 we allow an additional scalar Kit amplitude, the " K , as a 
Breit-Wigner with mass dependent width. Gaussian form factors are used for i^'Q(1430) and 
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K. For all Ktt resonances with non-zero spin, the radii in the Blatt-Weisskopf [16(] form factors 
are fixed, td = 5 GeV~^, = 1.5 GeV~^. Fit fractions obtained in our fit are compared 
with E791 in Table [D We get a poor fit quality xVi^=448/388, P{x^,u)=2%; the NR fit 
fraction is ~10%; and the kit fit fraction is dominant, ~30%. For the ii'g(1430) resonance 
we measure ^^.(1430) = (1461 ± 3) MeV/c^ and ri^.(i43o) = (169 ± 5) MeV/c^, which are 
consistent with E791 results, but inconsistent with current PDG |l5| values. Similar behavior 
is reported in a recent preprint from the FOCUS collaboration [l7|. We also measure the 
K, resonance parameters, = (805 ± 11) MeV/c^, and = (453 ± 21) MeV/c^. In these 

TABLE I: Comparison of fit fractions in percent between this work and E791 results for an isobar 
and a QMIPWA model. This work is preliminary and only statistical errors are shown. 





Model C 


QMIPWA 


Mode 


E791 f5] 


CLEO-c 


E791 [6] 


CLEO-c 


NR 


13.0±5.8±4.4 


10.4±1.3 


see S wave 


see S wave 


K*(892)7r+ 


12.3±1.0±0.9 


11.2±1.4 


11.9±0.2±2.0 


10.0±0.3 


Ko(1430)^+ 


12.5±1.4±0.5 


10.5±1.3 


see S wave 


11.4±3.6 


K2(1430)7r+ 


0.5±0.1±0.2 


0.40±0.04 


0.2±0.1±0.1 


0.476±0.014 


K*(1680)7r+ 


2.5±0.7±0.3 


1.36±0.16 


1.2±0.6±1.2 


2.52±0.08 




47.8±12.1±5.3 


31.2±3.6 


see S wave 


see S wave 


Total S wave 


73±15 


52±4 


78.6±1.4±1.8 


67.4±1.3 


x'^/u, Prob.(%) 


46/63, 94% 


448/388, 2% 


277/277, 47.8% 


368/346, 19.5% 



three models our fits well reproduce behavior reported by the E791 [5l] experiment. We get 
results for phases, fit fractions, and resonance parameters which are statistically consistent 
with E791. However the poor 2%-probability of the fit motivates us to explore an alternative 
model of the decay amplitude. 

In an approach motivated by Ref. [oj we set the form factor for any S-wave contribution 
to unity. Following prescriptions from Ref. [13] we replaced the Breit-Wigner function for 
K by the complex pole approximation WK(m) = [m^ — m^]~^, which is equivalent to the 
Breit-Wigner with a fixed width. We find pole parameters, S^m^ = (651 ± 17) MeV/c^ and 
Q^m^ = (—229 ± 22) MeV/c^, which are consistent with the values obtained in Refs. 
and pill]. We find no evidence of production or fit improvement by adding a i^*(1410)7r"'" 
intermediate state. We also tested the hypothesis j§| that the i^Q(1430) shape in this decay 
might be distorted by the vicinity of the Krf threshold. We replace the Breit-Wigner 
function for _ft'Q(1430) by the Flatte function [l8| : 

>Vx*(1430)("^) = ["^^.(1430) -m'^ - i{gK*K^pK^{m) + g\*KripKr,{m) + ^X*Xr,'Pi^'7' 

where g'^ and p(m) = 2P/m are the coupling constants and the phase space factors, re- 
spectively. We get the values of Flatte parameters, ruK'^mo) = (1468.3 ± 4.2) MeV/c^, 
gK'Kn = (543±11) MeV/c^, gK*Kr]' = (268±86) MeV/c^, and gx^Kri statistically consistent 
with zero, and the fit quality is not meaningfully improved. 

In all of our fits that contain only Kn resonant contributions, we observe that the p.d.f. 
significantly deviates the data in the range 1.3 < m^(7r"'"7r"'") < 1.8 (GeV/c^)^, as shown in 
Fig [31 As a solution to this problem we include an isospin-two 7r+7r+ S-wave contribution 
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FIG. 3: The m2(7r+7r+) projec- FIG. 4: The Dahtz plot projections in QMIPWA. 

tion in Model C. 
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FIG. 5: The Kt\ S wave in the isobar model (IM) and QMIPWA. 

in unitary form from Ref. 3] which accounts for the mass-dependent phase measured in 
scattering experiments 23]. With its inclusion, we find up to 20% improvement of the 
fit probability if we allow the inelasticity parameter to drop from unity to a much lower 
value above the p — p threshold starting from m^(7r+7r+) = 1.245 GeV/c^, using the smooth 
threshold function, as in Ref. [21]. 

The framework of the isobar model is good for a first approximation to data and compar- 
ison with other experiments. However, this model violates basic principles, such as unitarity 
and analyticity. It is useful to provide a model independent measurement of the partial 
waves. E791 in Ref. presents the first quasi model-independent partial wave analysis 
(QMIPWA) of the — »• K~'j:^i^^ decay. We apply the same technique with minor modi- 
fications. 

To describe the K-n amplitude and phase in a model independent way we use binned 
complex amplitudes. The entire range of vn^iKv:) from 0.4 to 3 (GeV/c^)^ is split into 
26 uniform bins, with width 0.1 (GeV/c^)^. Our fit has 26 amplitudes and 26 phases as 
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FIG. 6: The K-k P wave in the isobar model 
(IM) and QMIPWA. 



FIG. 7: The K-k D wave in the isobar model 
(IM) and QMIPWA. 
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FIG. 8: The 1=2 7r+7r+ S wave in the isobar 
model and QMIPWA. 



FIG. 9: The Kii S wave in the isobar model 
and QMIPWA with additional 1=2 7r+7r+ S- 
wave amplitude. 



fit parameters. Similarly, the tttt amplitude and phase is defined by 18 bins in of m^(7r7r) 
from 0.1 to 1.9 (GeV/c^)^. Wc use linear interpolation between bin centers. To alleviate 
the problem of coarse interpolation we use Brcit-Wigners for parameterization of "narrow" 
resonance structures. Thus, our S wave is made up of the Flatte function for the i^rQ(1430) 
and a binned S-wave amplitude, which replaces the k and non-resonant contributions in 
the isobar model. In contrast to E791 we do not use form factors for S-wave contributions. 
Our P wave is made up of the Breit-Wigner function for the ir*(892) and a binned P 
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amplitude, which substitutes for the wide i^*(1680) Breit-Wigner function in the isobar 
model. Our D wave contains a binned D amplitude which substitutes for the J^gll^SO) 
Breit-Wigner function of the isobar model. The isospin-two vrvr S wave, in addition to the 
unitary amplitude, also has a binned amplitude. 

The parameters of the binned functions describing the S, P, and D waves are allowed to 
float one wave at a time. Other resonance parameters for the same wave are fixed to their 
optimal values from the isobar model. The resonance parameters of the other spin waves 
float, as in the isobar model. First we fit with a floating binned amplitude for the S wave, 
using the isobar model approximation for P and D waves. The results of this fit are shown in 
Fig. [5] and Tabled We find an almost constant S-wave binned amplitude in the complete Kt: 
phase space, which differs slightly from substituting isobar model components. The binned 
phase shows a smooth rise from —120° to —10° from minimum to maximum invariant mass. 
The total S-wave amplitude is distorted by the i^'o(1430) resonance, as shown by the solid 
curve. 

Next, we fix parameters of the binned S wave and allow the P or D waves to float. The 
result of these fits and their comparison with the isobar model are shown in Figs. [H] and [3, 
respectively. We find that in both cases the P and D binned amplitudes are consistent with 
the results of the isobar model. For the D wave amplitude the resolution is poor, because 
its fraction is well below 1%. 

In a fit with floating binned amplitudes for the isospin-two 7r+7r+ S wave we also find a 
consistent description with the analytical function, shown in Fig. [HI Adding the isospin-two 
amplitude in its analytical form the fit with the Kit S wave slightly changes the results for 
the binned Ktx S wave, as shown in Fig. O 

In systematic cross checks we test the stability of our results by varying the event selection, 
efficiency and background parameterization, variation of model parameters and resonant 
components, etc. We expect that systematic uncertainties will be of the same size as the 
statistical errors. The results of this CLEO analysis are preliminary, and we plan to publish 
them after completion of systematic cross checks. 

In summary, we present the Dalitz plot analysis of the — >■ K~t[^tx^ decay using cur- 
rently available CLEO-c data. Using a simple isobar model we find that our data are consis- 
tent with results obtained in E791 experiment [sj. Inclusion of the isospin-two 7r+7r+ S wave 
improves consistency of this model with data. Finally, we apply quasi model-independent 
partial wave analysis [6^] and measure the amplitude and phase of the Ktt and 7r"^7r+ S waves 
in the range of invariant masses from the threshold to the maximum in this decay. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the A. P. Sloan Foundation, 
the National Science Foundation, the U.S. Department of Energy, and the Natural Sciences 
and Engineering Research Council of Canada. 
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